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We present measurements of the branching fraction and time-dependent CP asymmetries in 
B° -> J/ipir° decays based on (231.8 ± 2.6) x 10 6 T(4S) -* BB decays collected with the BABAR 
detector at the SLAC PEP-II asymmetric-energy B factory. We obtain a branching fraction B(B° — > 
J/^tt ) = (1.94 ± 0.22 (stat) ± 0.17 (syst))x 10 -5 . We also measure the CP asymmetry parameters 
C = -0.21 ± 0.26 (stat) ± 0.06 (syst) and S = -0.68 ± 0.30 (stat) ± 0.04 (syst). 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Charge conjugation-parity (CP) violation in the B 
meson system has been established by the BABAR [1] 
and Belle [2] collaborations. The Standard Model (SM) 
of electroweak interactions describes CP violation as a 
consequence of a complex phase in the three-generation 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix [3] . Measurements of CP asymmetries in the proper- 
time distribution of neutral B decays to CP eigen- 
states containing a charmonium and K° meson pro- 
vide a precise measurement of sin 2/3 [4], where (3 is 

^gE-KdKb/^td^tb] and the V U are CKM matrix el- 
ements. 

The decay B° -> J/tpn is a CP-even Cabibbo- 
suppressed b — * ccd transition whose tree amplitude has 
the same weak phase as the b — > ccs modes e.g. the CP- 
odd decay B° — > J/ipK®. The b — > ccd penguin ampli- 
tude has a different weak phase than the tree amplitude. 
The tree and penguin amplitudes expected to dominate 
this decay are shown in Figure 1. 



B° 
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FIG. 1: Feynman diagrams of the color suppressed tree (top) 
and gluonic penguin (bottom) amplitudes contributing to the 
B° -v J/077 decay. 



If there is a significant penguin amplitude in B° — > 
J/ipTT°, then one will measure values of the CP asymme- 
try coefficients S and C that are different from — sin 2(3 



and 0, respectively [5]. The coefficient S denoting the in- 
terference between mixing and decay, and the direct CP 
asymmetry coefficient C are defined as: 

22mA f-|A| 2 

S ^TTW and Ce IW' (1) 

where A is a complex parameter that depends on both 
the B°-B° oscillation amplitude and the amplitudes de- 
scribing B° and B° decays to the J/ipir final state. An 
additional motivation for measuring S and C from B° — > 
J/0 7T° is that they can provide a model-independent con- 
straint on the penguin dilution within B° —> J/0 Kg [6] . 

In this publication, we present an update of previ- 
ous BABAR branching fraction and time-dependent CP 
violating asymmetry measurements of the decay B° — » 
J/07T [7, 8], which had been performed using 20.7 fb _1 
and 81.1 fb _1 of integrated luminosity, respectively. Belle 
has also studied this mode and has published a branching 
fraction and later a time-dependent CP violating asym- 
metry result using 29.4 fb _1 and 140.0 fb -1 of integrated 
luminosity, respectively [9, 10]. 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric e + e~ stor- 
age ring. This represents a total integrated luminosity of 
210.6 fb _1 collected on or just below the T(4S) resonance 
(on-peak), corresponding to a sample of 231.8 ± 2.6 mil- 
lion BB pairs. An additional 21.6 fb" 1 of data, collected 
at approximately 40 MeV below the T(4S) resonance, is 
used to study background from e + e~ — > qq (q = u, d, s, c) 
continuum events. 

The BABAR detector is described in detail else- 
where [11]. Surrounding the interaction point is a 5 layer 
double-sided silicon vertex tracker (SVT) which measures 
the impact parameters of charged particle tracks in both 
the plane transverse to, and along the beam direction. A 
40 layer drift chamber (DCH) surrounds the SVT and 
provides measurements of the transverse momenta for 
charged particles. Both the SVT and the DCH oper- 
ate in the magnetic field of a 1.5 T solenoid. Charged 
hadron identification is achieved through measurements 
of particle energy loss (dE/dx) in the tracking system and 
the Cerenkov angle obtained from a detector of internally 
reflected Cerenkov light (DIRC). This is surrounded by 
a segmented CsI(Tl) electromagnetic calorimeter (EMC) 
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which is used to provide photon detection and electron 
identification, and is used to reconstruct neutral hadrons. 
Finally, the instrumented flux return (IFR) of the magnet 
allows discrimination of muons from pions. 

We reconstruct B° — > J/ipn° decays in BB candi- 
date events from combinations of J/ip — > C + Cr (I = 
e, /i) and ir° — > 77 candidates. A detailed descrip- 
tion of the charged particle reconstruction and identi- 
fication can be found elsewhere [7]. For the J/ip — > e + e~ 
{J/ip — > fi + /J.~) channel, the invariant mass of the lep- 
ton pair is required to be between 3.06 and 3.12GeV/c 2 
(3.07 and 3.13 GeV/c 2 ). Each lepton candidate must also 
be consistent with the electron (muon) hypothesis. We 
form 7T° — > 77 candidates from clusters in the EMC with 
an invariant mass, m 77 satisfying 100 < to 77 < 160 
McV/c 2 . These clusters are required to be isolated from 
any charged tracks, carry a minimum energy of 30McV, 
and have a lateral energy distribution consistent with 
that of a photon. Each ir° candidate is required to have 
a minimum energy of 200 MeV and is constrained to the 
nominal mass [12]. Finally the B° — > J/ip it candidates 
(B rec ) are constrained to originate from the e + e~ inter- 
action point using a geometric fit. 

We use two kinematic variables, Toes and AE, in order 
to isolate the signal: m ES = y/ (E^ cam ) 2 - (p* B ) 2 is the 
beam-energy substituted mass and AE = E* B — E^ cam is 
the difference between the -B-candidate energy and the 
beam energy. E^ cam and p* B {E* B ) are the beam energy 
and .B-candidate momentum (energy) in the center-of- 
mass (CM) frame. We require Toes > 5.2GeV/c 2 and 
I AS I < 0.3 GeV. 

A significant source of background is due to e + e~ — > 
qq (q = u,d,s,c) continuum events. We combine sev- 
eral kinematic and topological variables into a Fisher 
discriminant (J 7 ) [13] to provide additional separation 
between signal and continuum. The three variables 
L , L 2 and cos(0h) are inputs to T. L and L 2 are 
the zeroth- and second-order Legendre polynomial mo- 
ments; L = £ 2 |p*| and L 2 = £. | p * |/2 (3 cos 2 6 t - 1), 
where p* are the CM momenta of the tracks and neutral 
calorimeter clusters that are not associated with the sig- 
nal candidate. The 9i are the angles between p* and the 
thrust axis of the signal candidate and 9h is the angle be- 
tween the positively charged lepton and the B candidate 
in the J/ip rest frame. 

We use multivariate algorithms to identify signatures 
of B decays that determine (tag) the flavor of the de- 
cay of the other B in the event (B tag ) to be either a B° 
or B°. The flavor tagging algorithm used is described 
in more detail elsewhere [14]. In brief, we define seven 
mutually exclusive tagging categories. These are (in or- 
der of decreasing signal purity) Lepton, KaonI, Kaonll, 
Kaon-Pion, Pion, Other, and No- Tag. The total effective 
tagging efficiency of this algorithm is (30.5 ± 0.4)%. 

The decay rate / + (/_ ) of neutral decays to a CP eigen- 



state, when B tag is a B° (B°), is: 

e -\At\/r B o 

f±(At) = — [l±5 , sin(Am d At)=FCcos(Am d At)], 

4t b o 

(2) 

where At is the difference between the proper decay times 
of the B rec and B tag mesons, t b o — 1.536 ± 0.014 ps is 
the B° lifetime and Am d = 0.502 ± 0.007 ps" 1 is the B°- 
B° oscillation frequency [12]. The decay width difference 
between the B° mass eigenstates is assumed to be zero. 

The time interval At is calculated from the measured 
separation Az between the decay vertices of B rec and 
B tag along the collision axis (z). The vertex of B rec is 
reconstructed from the lepton tracks that come from the 
J/ip; the vertex of B tag is constructed from the remain- 
ing tracks in the event that do not belong to B rec , with 
constraints from the beam spot location and the B rec 
momentum. We accept events with \At\ < 20 ps whose 
uncertainty are less than 2.5 ps. 

After all of the selection criteria mentioned above have 
been applied, the average number of candidates per event 
is approximately 1.1, indicating some events still have 
multiple candidates. In these events, we randomly choose 
one candidate to be used in the fit. This selection is 
unbiased. Overall, the true signal candidate is correctly 
identified 91.7% of the time. After this step, the signal 
efficiency is 22.0% and a total of 1318 on-peak events are 
selected. 

In addition to signal and continuum background 
events, there are also £?_B-associated backgrounds present 
in the data. We divide the B backgrounds into the follow- 
ing types: (i) B° — > J/ipK®, where — > 7r°7r° (ii) inclu- 
sive neutral B meson decays, and (iii) inclusive charged 
B meson decays. When normalized to the integrated 
luminosity, Monte Carlo (MC) studies predict 153 ± 9, 
68 ± 14 and 314 ± 63 events of these background types, 
respectively. The inclusive neutral B meson decays ex- 
clude signal and B° — > J/ip Kg events. The inclusive 
B decay backgrounds are dominated by contributions 
from B — > J/ipX (inclusive charmonium final states). 
In particular the inclusive charged B meson decay back- 
grounds are dominated by B ± — > J/tpp + decays. The 
B° —> J/ip K® background was studied separately since 
its CP asymmetries are precisely measured. 

The signal yield, S and C are simultaneously extracted 
from an unbinned maximum-likelihood (ML) fit to the 
B candidate sample, where the discriminating variables 
used in the fit are Toes, AE, T and At. The continuum 
yield for the seven mutually-exclusive tagging categories, 
is also allowed to vary in the ML fit. 

The probability density function (PDF) for signal Toes 
distribution takes the form of a Gaussian with a low side 
exponential tail [15]. We parameterize the toes distribu- 
tion for continuum and neutral inclusive B background 
with an Argus phase space distribution [16]. As there 
are significant correlations between toes and AE for the 
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charged inclusive B and the B° — > J/ipK® backgrounds, 
we parameterize these variables with two-dimensional 
non-parametric PDFs. The AE distribution for signal 
events is modeled by a Gaussian with an exponential tail 
on the negative side to account for energy leakage in the 
EMC, plus a polynomial contribution. The AE distribu- 
tions for the continuum and the neutral inclusive B back- 
ground are described by second and third-order polyno- 
mials respectively. The T distributions for the signal and 
the backgrounds are described by bifurcated Gaussians 
with different widths above and below the peak value. 

The signal decay rate distribution of Equation 2 is 
modified to account for dilution coming from incorrectly 
assigning the flavor of B tag and is convolved with a triple 
Gaussian resolution function, whose core width is about 
f.l ps [17]. The decay rate distribution for B back- 
grounds is similar to that for signal. The inclusive B 
backgrounds are assigned an effective lifetime instead 
of their respective B lifetimes to account for their mis- 
reconstruction. This effective lifetime is determined from 
MC simulated data. The potential CP asymmetry of the 
inclusive B background is evaluated by allowing the pa- 
rameters of S and C for this background to vary. The 
decay rate distribution for B° — > J/tpK® is the same as 
that for signal and reflects the known level of CP viola- 
tion in that decay. The continuum background is mod- 
eled with a prompt lifetime component convolved with a 
triple Gaussian resolution function. The core Gaussian 
parameters and fractions are allowed to vary in the ML 
fit. The other two Gaussians have means fixed to zero, 
and widths of 0.85 ps and 8.0 ps, respectively. 

The results from the ML fit are 109 ± 12 (stat) sig- 
nal events, with S = -0.68 ± 0.30 (stat) and C= -0.21 
± 0.26 (stat). The fit yields the following numbers of 
continuum events: Ni epto „ = 17 ± 5, ^Kaoni — 38 ± 

8, Nxaonll — 101 ± 12, ^KaonPion = 102 ± 12, Np ion 

= 115 ± 12, Nother = 94 ± 11 and N NoTag = 227 ± 
17. Figure 2 shows the distributions of toes, AE, and 
T for the data. In these plots the signal has been en- 
hanced by selecting \AE\ < 0.1 GeV for the Toes plot, 
toes > 5.275 GeV/c 2 for the AE plot and by applying 
both of these criteria for the T plot. After applying 
these requirements to the signal (background) samples 
that are used in the fit, they are reduced to a relative 
size of 83.1% (24.3%), 85.0% (21.1%) and 73.1% (2.8%) 
for the toeSi AE, and T distributions respectively. 

Figure 3 shows the At distribution for signal B° 
and B° tagged events. The signal has been enhanced 
using the same toes and AE cuts as for Figure 2. 
The time-dependent decay rate asymmetry [N(At) — 
N(At)]/[N(At) +N(At)} is also shown, where N (N) 
is the decay rate for B° (B°) tagged events and the de- 
cay rate takes the form of Equation 2. 

Table I summarizes the systematic uncertainties on the 
signal yield, S and C. These include the uncertainty 
due to the PDF parameterization (including the resolu- 
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FIG. 2: Signal enhanced distributions of rriES (top), AE (cen- 
ter) and T (bottom) for the data (points). The solid line 
represents the total likelihood, the dashed line is the sum of 
the backgrounds and the dotted line is the signal. The un- 
dulations in the background model are the result of limited 
MC statistics available for defining the two-dimensional non- 
parametric PDFs. 



tion function), evaluated by varying the signal and the 
background PDF parameters within uncertainties of their 
nominal values, the effect of SVT mis-alignment; the un- 
certainties associated with the Lorentz boost, the z-scale 
of the tracking system, and the event-by-event beam spot 
position. 

The uncertainty coming from the fit bias is estimated 
by performing ensembles of mock experiments using sig- 
nal MC which is generated using the GEANT4-based [18] 
BABAR MC simulation, embedded into MC samples of 
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FIG. 3: The At distribution for a sample of signal enhanced 
events tagged as B° (top) and B° (middle). The dotted lines 
are the sum of backgrounds and the solid lines are the sum 
of signal and backgrounds. The time-dependent CP asymme- 
try (see text) is also shown (bottom), where the curve is the 
measured asymmetry. 



TABLE I: Contributions to the systematic errors on the sig- 
nal yield, S and C, where the signal yield errors are given in 
numbers of events. The total systematic uncertainty is the 
quadratic sum of the individual contributions listed. Addi- 
tional systematic uncertainties that are applied only to the 
branching fraction are discussed in the text. 



Contribution 


Signal yield 


S 


C 


PDF parameterization 


+3.21 
-2.88 


±0.013 


±0.012 


SVT mis-alignment 




±0.002 


±0.002 


Boost and z-scale 


+0.08 
-0.16 


±0.004 


±0.001 


Beam spot position 




±0.007 


±0.002 


Fit bias 


±3.00 


±0.026 


±0.016 


Inclusive B background yields 


±3.52 


±0.003 


±0.020 


rriEs-AE correlations 


±2.92 


±0.020 


±0.002 


CP content of B background 


+0.13 
-0.11 


±0.012 


±0.049 


CP background lifetime 


±0.67 


±0.010 


±0.010 


Tagging efficiency asymmetry 


±0.02 


±0.000 


±0.020 


Tag-side interference 




±0.004 


±0.014 


Fisher data/MC comparison 


±0.70 


±0.004 


±0.004 


Total 


+6.42 
-6.26 


±0.040 


±0.063 



background generated from the likelihood. The deviation 
from input values is added in quadrature to the error on 
the deviation in order to obtain a conservative fit bias un- 
certainty. Most, but not all of the inclusive charmonium 
final states which dominate the inclusive B background, 
are precisely known from previous measurements. Their 
yields are then fixed in the fit. 



As a crosscheck, the yields for inclusive B backgrounds 
that are not well known are allowed to vary. The devi- 
ation from the nominal result is taken as a systematic 
uncertainty. We include an additional systematic uncer- 
tainty to account for neglecting the small correlation be- 
tween tties and AE in signal and neutral inclusive B 
background events. 

In order to evaluate the uncertainty coming from CP 
violation in the B background, we have allowed the S 
and C parameters to vary in a fit for the neutral inclu- 
sive B background, and have separately allowed the C 
parameter to vary in a fit for the charged inclusive B 
background. The deviations of the fitted values of the 
signal S and C from the nominal fit results are assigned 
as systematic errors. The uncertainty from CP violation 
in B° — > J/ip K® is determined by varying S and C within 
current experimental limits [14]. 

The inclusive B background uses an effective lifetime 
in the nominal fit and we replace this with the world- 
average B lifetime [12] to evaluate the systematic error 
due to the CP background lifetime. There is also a small 
asymmetry in the tagging efficiency between B° and B a 
tagged events, for which a systematic uncertainty is eval- 
uated. We study the possible interference between the 
suppressed b — > ucd amplitude with the favored b — > cud 
amplitude for some tag-side B decays [19]. The difference 
in the distribution of T between data and MC is evalu- 
ated with a large sample of B — > D* p decays. There are 
additional systematic uncertainties that contribute only 
to the branching fraction. These come from uncertain- 
ties for charged particle identification (5.2%), ir° meson 
reconstruction efficiency (3%), the J/ip —> £ + £~ branch- 
ing fractions (2.4%), the tracking efficiency (1.2%) and 
the number of B meson pairs (1.1%). The systematic 
error contribution from MC statistics is negligible. The 
109 ±12 signal events correspond to a branching fraction 
of 

B(B° -> J/V>7r°) = (1.94±0.22(stat)±0.17(syst)) x 1CT 5 . 

We determine the CP asymmetry parameters to be 
C = -0.21±0.26(stat)±0.06(syst), 
S = -0.68 ± 0.30(stat) ± 0.04(syst), 

where the correlation between S and C is 8.3%. The 
value of S is consistent with SM expectations for a tree- 
dominated b — ► ccd transition of S = — sin 2/3 and C = 
0. All results presented here are consistent with previous 
measurements from the B Factories [7-10]. 
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